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By using an approach based on the Self-Consistent Local Random Phase Approximation (SC- 

LRPA) we calculate the magnetic excitation spectrum in diluted/disordered ferromagnetic systems. 
In previous studies, the SC-LRPA was shown to be reliable for the determination of the Curie 
temperature. In this letter, we now demonstrate its accuracy and efficiency for the determination of 
the magnetic excitation spectrum. For that purpose, we calculate the spin excitations of the widely 
studied diluted III-V magnetic semiconductor Gai-^Mn^As. Instead of the expected broadening of 
the excitations due to disorder, it is shown that magnons exist only in a very restricted region of the 
Brillouin zone. We calculate without adjusting parameters, the spin stiffness in optimally annealed 
systems as a function of x and compare it to recent measurements. We find a very good agreement 
for well annealed samples and provide a plausible explanation for the very small values measured in 
as grown ones. We hope that this study will stimulate new studies as Ineslatic Neutron Scattering 
measurements for example. 



The understanding of the influence and effects of dis- 
order on transport and magnetic properties in magnetic 
materials is a key issue from both a fundamental point 
of view or in the prospect of possible technological ap- 
plications. This requires the necessity to develop new 
methods and tools which are able to treat the disorder 
in a reliable manner. The word disorder is very gen- 
eral and includes for example, (i) the dilution of magnetic 
atoms in a non magnetic matrix, (ii) substitution of a 
non magnetic atom by another one which has a different 
radius or/and valency as in manganites Ai-^B^MnOa 
where A=La,Pr and B=Sr,Ca,Ba for example, (iii) in- 
trinsic defects which may appear during the growth of 
the material, for instance vacancies in compound 
as ZrC>2, HfG-2 or TiC>2. By browsing the literature, one 
can find several examples which underline the importance 
and crucial role of the disorder. For example, in mangan- 
ites one can observe after substitution, metal-insulator 
phase transition, formation of nanoscale inhomogeneities, 
or anomalies in the magnetic excitation spectrum, see for 
example 0, 0, S 0] • We can also mention the high Curie 
temperature reported in materials which, a priori, do not 
contain any or at least a sufficient amount of magnetic 
impurities Q. The ferromagnetism is attributed to 
defects which can in the case of HfC>2 appear close to the 
interface between the material and the substrate. An- 
other example is the spin wave ferromagnetic excitation 
observed after vacancies are introduced in the cobaltate 
Nai_ x D x Co02 III (the symbol □ corresponds to vacan- 
cies), thus the substitution of Na by a vacancy is more 
than just a way to vary the carrier concentration. In 
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diluted magnetic semiconductors as GaMnAs it is seen 
that the magnetic moment distribution and Curie tem- 
perature are also very sensitive to the sample history and 
thus to the presence of compensating defects as Mn in- 
terstitials and As anti-sites [8|. 

In this manuscript we focus our attention on the nature 
of the magnetic excitation spectrum in dilute ferromag- 
netic materials. We will show that the proper and simul- 
taneous treatment of both disorder and thermal fluctua- 
tions leads to unusual excitation spectrum in comparison 
to what is usually observed in non dilute systems. The 
paper is organized as follows. In the first part, we present 
a summary of the formalism based on the Self-Consistent 
Local Random Phase Approximation theory (SC-LRPA) . 
In the second part, using the exchange couplings calcu- 
lated within first principle approach (no adjustable pa- 
rameters), we calculate the spectral function, excitation 
spectrum and spin stiffness as a function of the impurity 
concentration in optimally annealed Gai-^Mn^As. The 
term optimally annealed means that the concentration 
of compensating defects is negligible. In the last part 
we compare our results to available experimental data 
for the spin stiffness and explain the large difference ob- 
served between as grown and annealed samples. 

The Hamiltonian which describes N imp interacting 
quantum (or classical) spins randomly distributed on a 
given lattice of N sites is the dilute Heisenberg model, 

u ;>>..s,.s, a) 

The couplings Jy are very general (short or long 
range), they are assumed to be given and can eventu- 
ally depend on the impurity concentration x — ^ v . 
We will come back to this point in what follows. For a 
given configuration of disorder e.g. the coordinates of the 
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magnetic impurities, we define the retarded spin Green's 
function G£(w) =' /+~ G%{t)e wt dt where G£(t) = 
-^(^([^(i),^^)]) which describes the transverse 
spin fluctuations, where the subscript "c" corresponds 
to the configuration index. After performing the Local 
Random Phase Approximation decoupling of the higher 
order Green's functions which appear in the equation of 
motion oiGUu) HGjJ, 

we find, 



where A£(q) = ^_ £y A i (<|*«' c >(*S' c |i)e <q(r '- r ^, 
we have introduced temperature dependent local param- 



eter A ,, 



In the following we will calculate di- 
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rectly A(q,uj). However, It is also interesting to define 
the moments associated to it, 
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m n(q) = / w"A(q, w)dw 

J — oo 
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where H£ ff , G c and D are Ni mp x Ni mp matrices. The 
index i runs over only the sites occupied by a local- 
ized spin. The effective Hamiltonian matrix elements is 
(H e c ff )„ = -{St)J l3 +5 l3 Zi(Sf)Jl, and D tj = 2<5f>%. 
The local magnetization (5f) has to be calculated self- 
consistently at each temperature. Note also that the con- 
dition J2j(K-lf[)ij = 0, implies that zero is eigenvalue trum depends on the renormalized temperature depen- 
dent couplings by J?f(T) — JijXjXi. Since A, are fluc- 



We now make some general remarks. It is 
straightforward to see that the first moment can 
be rewritten mi(q) = J(0) — J(q) where J(q) = 
((S z )) (^EyA.-Ai^-e^-^))^ This expression is 
interesting and shows that the magnon excitation spec- 



of H£ ff , the SC-LRPA treatment is consistent with the 
Goldstone theorem. Note that although the matrix is 
non Hermitian, in the ferromagnetic phase the spectrum 
is real and positive at each temperature. If a negative 
eigenvalue would appear, this would indicate an insta- 
bility of the ferromagnetic phase, as in the case where 
frustration would be present. In the following, we will 
discuss the case of GaMnAs for which all the couplings 
are ferromagnetic, thus such instabilities will not occur. 
Because H^ ff is non Hermitian (real non symetric), we 
now precise how the GF is expressed. H£ ff , has the prop- 
erty to be bi-orthogonal 14]. Thus, we define right and 
left eigenvectors of iJ c ff denoted respectively \^a' c ) an< ^ 
|^^' c ), both are associated to the same eigenvalues u>%: 
*5 C ) = uim*' c ) and *h!Lf l^' c > = u£l*£' c ). In 



tuating from site to site, the previous expression in- 
dicates that the excitation spectrum will have a more 
complex temperature dependence than that obtained in 
the non dilute case (x — 1 and A, = 1) for which 
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,{T = 0). In the limit of q -> 0, we 
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effl . a ) = wSl*£ c > and *H^|*^ C > = co c a \^ c ). 
general, two eigenvectors belonging to the same set L or 
R are not orthogonal to each other, but when they are of 
different type they fulfill the relation ' c |*«; c ) = 5 a>a >. 
After inserting the L,R eigenvectors in eq.©, the re- 
tarded Green's function Gfj (w) can be rewritten, 



immediately find that mi(q) « Diq 2 where the effective 
stiffness (in the following we will clarify the use of "effec- 
tive") is D X = ((5*)) Y,ij A , A,. /,,;,-, - r,) 2 ) c . 

Note that, at low temperature where all spins are po- 
larized, this expression can be simplified and does not 
depend on the nature of the eigenstates of the effective 
Hamiltonian. If we assume that all couplings are ferro- 
magnetic (as in GaMnAs) then we get, 
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where Z?i 
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where the double sum 
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Although the system is non translation invariant 
we define the Fourier transform by G(q, to) 
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r ^G\AuA) c . The notation („..) c de- 



notes the average over the disorder configurations. We 
now define the dynamical spectral function A(q, uS) — 
^^ImG^u) where <(5*» = ^ £,(5?) is the av- 
erage magnetization over the impurity sites. The spectral 
function A(q, uj) is the physical quantity which provides 
the direct access to the magnetic excitation spectrum, it 
is accessible by Inelastic Neutron Scattering experiment. 
It reads, 
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runs now over all sites. This expression indicates that the 
spin stiffness extracted from the first moment calculation 
reduces only to the virtual crystal expression. The conse- 
quence of this is the following. First D\(x) overestimates 
the real spin stiffness. Additionally, at sufficiently low 
dilution, one expects that below the percolation thresh- 
old both Tc and the spin stiffness should vanish. This 
expression is finite below the percolation threshold and 
does not know about it. Thus D\ does not correspond 
to the real spin stiffness. In the following we will explain 
how to evaluate it properly. 

The other interesting quantity is the effective linewidth 
7(q) = v fi2(q) — mf(q), it measures the broadening of 
the magnetic excitations due to the disorder (dilution) 
and thermal fluctuations. In the long wave length limit, 
we can show that 7(q) = Cq. The linewidth is larger 
than the excitation energy. Similar results were reported 
in the study of the disordered double exchange model 
(nearest neighbor coupling) relevant for manganite com- 
pounds (x = 1) [151 ]. Then, one would naively be tempted 
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FIG. 1: (Color online) Spectral function for Gai-^Mn^As as 
a function of the energy w (in Kelvin) in the (1,0,0) direc- 
tion for different values of q x , with x = 0.03. In the inset, 
are plotted the first moment, the effective linewidth and the 
position of the peak in A(q,w) as a function of q x /n(in units 
of ag 1 ) and the y-axis is in Kelvin. 



to conclude that there is no well defined magnetic excita- 
tions. As it will be seen in the following, even though the 
magnon excitations are incoherent, they appear as well 
defined peak in A(q, u>). In other words, the energy and 
width of the excitations calculated with the moments of 
the spectral function do not correspond to the peak posi- 
tion and linewidth of the excitations directly calculated 
with the full spectral function. 

In the following section and to illustrate our theory we 
calculate the magnetic excitation spectrum in the III-V 
dilute magnetic semiconductor Gai-^Mn^As as a func- 
tion of the Mn 2+ concentration and at T = K. In 
this case, A; = 1, the matrix H e g is real symetric and 
thus L and R eigenvectors are identical. We present a 
detailed analysis of the excitation spectrum as a func- 
tion of the Mn 2+ concentration. In addition, in order 
to allow a quantitative study with no fitting parameters 
and a direct comparison with the experimental measure- 
ments, the magnetic exchange coupling we use were cal- 
culated from first principle Tight Binding Linear Muffin 
Tin Orbital approach. The same couplings were used suc- 
cessfully for example in ref . [n| to calculate the Curie 
temperature. We stress that the used exchange couplings 
extend over 62 shells. 

In Figfl]we have plotted A(q,u>) for Gai-zMn^As as 
a function of energy for different values of the momen- 
tum q in the (1,0,0) direction. We remind that GaAs 
has a zinc blend structure, thus the magnetic impurities 
are distributed randomly on the fee sublattice of Ga. In 
addition, periodic boundary conditions have been used. 
The calculations were performed on a system which con- 
tains 4 x (32) 3 sites, the average over disorder was per- 
formed over 50 configurations. Note that for x = 0.03 
the system contains approximately Ni mp = 4000 impu- 



rities. In all cases, we have checked that the number of 
disorder configurations was sufficient to provide accurate 
results. Indeed, for example for the 4x (16) 3 system, it was 
found that the average over 50 configurations is already 
sufficient. Even if we increase the number of configura- 
tions A(q,ui) remains almost unaffected. 

In Fig[T]we can observe well defined excitations for only 
relatively small values of the momentum. As the mo- 
mentum q increases the peak becomes broader and the 
well defined excitation is then replaced by a very broad 
structure. In the inset we have plotted the first moment 
mi(q), the linewidth j(q) and the position of the peak 
which we denote as u>(q) as a function of the momentum 
q. It is clearly seen that at low q both mi(q) w D\q 2 
and u>(q) ~ Dq 2 . However the "effective stiffness" D\ 
appears to be larger than D, typically in this case we 
find Di/D « 2. This is in agreement with the previ- 
ous discussion that D\ overestimates the real spin stiff- 
ness, especially as we approach the regime of percola- 
tion. The main consequence of this is that the accurate 
determination of the spin stiffness should be done by di- 
rectly plotting the energy of the peak as a function of 
the momentum instead of using the first moment disper- 
sion. As mentioned in the previous section we indeed 
observe that the effective linewidth 7(17) is linear in q, 
but it does not correspond to the real linewidth of the 
peak observed in A(q,u>). The second moment is inap- 
propriate to evaluate the real linewidth of the magnetic 
excitation, this is reflected by the fact that the peak ob- 
served in FigJT] are non Gaussian-like and very asymmet- 
ric. In fact, asymmetric peaks in the magnetic excitation 
spectrum were observed in the diluted antifcrromagnet 
Mn x Zni_ x F2 ll|. It was observed that as we approach 
the percolation threshold x c — 0.25 the line shape of the 
low magnetic excitation peaks becomes strongly asym- 
metric and develops a long tail extending towards higher 
energy. This is clearly seen in Fig[TJ It was shown that 
the peaks can be fitted by a two component structure: 
a sharp Gaussian peak which describe the magnon exci- 
tation (its energy is q-dependent) and a broad localized 
mode (weakly q dependent). The observed increase of 
asymmetry of the peak as we increase the momentum 
corresponds to a cross over from propagating low energy 
excitations to localized excitations (fractons) (l2. 13|. 

In the next figure, Fig^l we have plotted the spectral 
function in (q, uS) plane. As in the previous figure the 
density of magnetic impurity is x = 0.03. In contrast to 
what is usually observed in weakly disordered magnetic 
systems as in manganites well defined excitations exist in 
the dilute case only in a restricted region of the Brilloum 
zone centered around the V point (q — (0,0,0)). This 
unusual feature shows that the nature of ferromagnetism 
in the diluted semiconductors is very different from that 
observed in non dilute materials. As seen in the previous 
figure the broadening of the excitation increases signif- 
icantly as we move away from the T point. We expect 
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FIG. 2: (Color online) Spectral function A(q, lu) in 
Gai-zMn^As in the (q,uj) plane (a: = 0.03). The energy 
axis (y axis) is in Kelvin 
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FIG. 4: (Color online) Spin stiffness D in meVA 2 as a function 
of the concentration x in optimally annealed samples and Tc 
in Kelvin in Gai-^Mn^As from ref. 0]. In the inset the ratio 
Tc/D (in A _2 )as a function of x is plotted. 
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FIG. 3: (Color online) Magnon energy u)(q) in meV as a func- 
tion of q 2 (in A~ 2 ) for different concentration of Mn 2+ im- 
purities (x = 0.02,0.03,0.04 and 0.05). The size L of the fee 
lattice varies from L = 16 to L = 32. 



that the momentum cut-off below which a well defined 
excitation exists should be related to the percolation cor- 
relation length £ p 



11] 



In FigJ3]we have plotted, for different concentration of 
Mn impurities, the spin wave energy to(q) at low q. To be 
more precise the energy of the first peak in A{q 1 u>) in the 
(1,0,0) (q = (j^-,0,0) is shown as a function of q 2 , the 
system size varies from L — 16 to L = 32. We clearly ob- 
serve that for small momentum q w(q) oc D(x)q 2 . Note 
that we have used for ao = 5. 65 A the value of the lattice 
spacing in GaAs.Note that, the error bars are included in 
the size of the symbols. If the average over disorder con- 



figuration was not sufficient one would expect that the 
curve u)(q) would exhibit some noise. In Fig HI we have 
plotted the spin stiffness D(x) and the Curie temperature 
(from ref.[lfj|) as a function of the Mn 2+ concentration 
x. We remind that there is no adjustable parameters in 
these calculations. We observe that the values of D(x) 
are relatively high for dilute ferromagnetic system. For 
example for x = 0.03 for which the Curie temperature is 
T c « 80 if the spin stiffness is D « 210 meVA 2 . Inter- 
estingly this value is close to what is usually measured in 
non dilute systems as in mang anites HE [H d ■ When 
approaching the percolation threshold from above the 
spin stiffness decreases very fastly, faster than the Curie 
temperature which can be fitted by Tc ~ A(x — Xc) 1 ^ 2 
[Io| . However we can not extract the critical power in 
a reliable manner, this should require more calculations 
close to the percolation threshold. In the inset we have 
also plotted the ratio T c /D, it is clearly seen that the 
Curie temperature is not proportional to the spin stiff- 
ness. 

It is now interesting to compare these results to ex- 
perimental data. In contrast to Curie temperature or 
resistivity measurements, there are only few experimen- 
tal studies concerning the magnetic excitation in diluted 
magnetic semiconductors. In our knowledge there are 
two papers in which the measurement of the spin stiffness 
in GaMnAs was done. In ref. the authors have mea- 
sured it by ferromagnetic resonance in as grown sample 
where x = 0.05. They have reported a value of approx- 
imately 60meV A 2 . In a recent study [2(j, done using 
samples of different concentrations and sample history, 
it was shown that annealed sample possess much higher 
values of the spin stiffness. Indeed, it was shown that the 
annealed sample (x = 0.03) which orders ferromagneti- 
cally below T c xp — 80K has a spin stiffness of the order 
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of D exp w WOmeVA 2 , whilst for as grown samples (x = 
0.05 and 0.06) D exp was only approximately AOmeVA 2 . 
We remind that our calculation are performed for opti- 
mally annealed samples. In FigfJ] for x = 0.03 we find 
Tg h = 88 K and D th (x = 0.03) = 2l0meVA 2 , both of 
these values agrees well with the experimental values. 
Because the concentration of Mn is not precisely known 
in experimental samples, a value of x = 2.5% would lead 
to Tq 1 = 75K and D th (x = 0.025) = 175meV"A 2 , this 
would significantly improve the agreement. Now the in- 
teresting question is: Why are the values of the spin stiff- 
ness is as grown sample so small? Let us assume that 
before annealing the main effect of compensation (for 
more details see ref.[l3|) is a reduction of the density of 
the magnetically active Mn (a; e ff). Because the measured 
T c « 65 - 80 K for x « 5% - 6% , from the theoretical 
curve of Tc we would conclude that x c s > 2.5%. Thus 
the expected stiffness should be larger than 150 meVA 2 . 
This is in disagreement with the much smaller experimen- 
tal values. We argue that the physical mechanism which 
appears to be responsible for such small values should 
be the inhomogeneous nature of the as grown sample in- 
stead of the reduction of the magnetically active localized 
spins. In ref. [l0|, we could explain the measured Tc in 
as grown sample by just estimating for each of them x e fi 
but with no need to take into account of the inhomoge- 
neous nature. Thus D is more sensitive to the presence 
of inhomogeneities than the Curie temperature [21| . 

To conclude, we have calculated the magnetic exci- 
tation spectrum in dilute ferromagnets within a theory 
based on the Self-consistent Local RPA. The theory is 
general and allows to treat disorder and thermal fluctu- 
ations in a reliable manner. As an illustration, we have 
calculated the spectrum in well annealed sample of GaM- 
nAs as a function of the impurity concentration. It is pre- 
dicted that well defined excitation can be observed only 
in a restricted region of the momentum space centered 
around the T point. A good agreement with experimen- 
tal values of the spin stiffness was obtained for the an- 
nealed sample with concentration x = 0.03. We have also 
discussed and provided an explanation to the very small 
value of the spin stiffness measured in as grown samples. 
We hope that this theoretical study will motivate new 
experimental measurements for example D as a function 



of the impurity concentration and as a function of the 
annealing treatment. 
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